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1 Introduction 



The EMC measurement of the (7i-structure function of the nucleon has given rise to 
renewed interest in the spin effects in high energy processes. In particular, many authors 
have been attracted by the chiral-odd spin-dependent distribution functions /ii.iP, H, 
Because of the chirahty, they can not be measured in the totally inclusive deep inelastic 
lepton- nucleon scattering (except as a quark-mass effect), but reveal themselves in polarized 
hadron-hadron collisions (Direct photon production, Drell-Yan process, etc.) and semi- 
inclusive polarized electron scatterings. Thus they are expected to open up a new window 
to explore hadron structure. 

The twist-2 parton distribution hi was first addressed by Ralston and Soper more than 
ten years ago and was recently named as the "transversity" distribution in ref. 0]. It 
has a simple parton model interpretation like /i and gi,^ and appears as a leading contribu- 
tion, for example, in the Drell-Yan process with the transversely polarized nucleon-nucleon 
collision . 

The other chiral-odd distribution hi is twist-3. Physically higher twist distributions 
represent complicated quark-gluon correlations in hadrons. It is difficult to isolate them by 
experiments because they are usually hidden by the leading twist-2 contributions. However, 
the twist-3 distribution hi is somewhat immune to this difficulty: hi reveals itself as a 
leading contribution to the longitudinal-transverse asymmetry in the Drell-Yan process |^ , 
although it appears with a factor 1 / ^/Q^. This is similar to the circumstance of g2 (the 
transverse spin structure function) whose contribution to the transversely polarized DIS 
becomes leading order but with the suppression factor 1 / y/Q^ [0 . Thus, the uniqueness of 

is twofold: it is a "measurable" higher-twist distribution and corresponds to chirality 
violating process. It is expected to provide new information about the hadron structure and 
the QCD dynamics beyond the conventional structure function data. 

In addition to more precise measurements of the familiar distribution functions f\ and 
gi, these "new" distribution functions hi, hL, and gj^ will be measured in the future col- 
lider experiments, such as HERMES, SMC ^ and RHIC. In view of this, it is especially 
important to develop theoretical study of these distribution functions as much as possible 
based on QCD. Among these efforts, the first step is the perturbative QCD prediction on 
the Q^-evolution of the distribution functions: Owing to the factorization property of hard 
processes, the Q^-evolution of the distribution functions can be predicted unambiguously 
in the framework of the renormalization group and the QCD perturbation theory. Its pre- 
diction is indispensable to extract physical information from the experimental data in the 
high-energy scale by comparing them with the prediction of hadron models at the low-energy 
scale. Furthermore, the comparison of the Q^-evolution itself between theory and experiment 
will provide a deeper test of QCD beyond the conventional twist-2 level. The Q^-evolution 
of the twist-2 distribution functions has been fully discussed since the first application of 
QCD to hard processes; for example, the Q^-evolution of hi was studied in p[ by employing 
the Altarelli-Parisi equation As for the higher twist ones, there has been several works 
on the Q^-evolution of gx by generalizing the Altarelli-Parisi equation to the higher twist 

■^We denote the twist-2 distributions corresponding to the structure functions Fi {F2) and gi by /i and 
gi following ref. Q . 

'^gr denotes the twist-3 distribution corresponding to the structure function (72- 
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distribution lITOl, ITT 



TIO, IT41IT5I1. But there 



and by the anomalous dimension calculation 
has been no discussion on the Q^-evolution of another important twist-3 distribution hi. 

In this paper we study the Q^-evolution of hi. We shall calculate the anomalous dimen- 
sion matrix for the twist-3 distribution hi based on the standard QCD perturbation theory. 
In general, the moments of a twist-3 distribution can be written in terms of the matrix ele- 
ments of a set of twist-3 operators involving explicitly the gluon field strength tensor, and the 
mixing among them occurs through renormalization, as was emphasized in ref. [|l^, Q in the 
context of the (72-structure function. However, the operator mixing occurs not only among 
these twist-3 operators but also with the other twist-3 operators which vanish by the naive 
use of the QCD equation of motion, {ip — mg)ip = 0, (referred to as the "equation-of-motion 
(EOM) operators" from now on). This is due to the fact that the naive equations of motion 
and thus the vanishing of these operators are not correct as an operator statement because 
of quantum effects and renormalization. The use of the equations of motion is allowed only 



when their matrix elements are taken with respect to a physical state [1T6| , |T^ . On the other 
hand, the renormalization of composite operators has to be carried out in terms of general 
Green functions which imbed these composite operators. Therefore the mixing involving the 
EOM operators is essential to perform renormalization of the higher twist operators consis- 
tently, which was recently pointed out by Kodaira, Yasui, and Uematsu |15] in the context of 
g2. We shall pay particular attention to this mixing and we will find that it certainly plays 
a role also for the present case of hi. 

The outline of this paper is as follows: In section 2, we shall first introduce the chiral- 
odd distributions hi and Hl, following the procedure of ref. [Q. We include this part to 
identify the manifestly interaction-dependent operators ("canonical basis" [16, 21 1) as well 
as the EOM operators relevant for the twist-3 part of hi- Readers familiar with ref. 
can skip this part by just noting the existence of the first term of eq. ( |2.10| ). Next we 
present the general procedure for the renormalization of the twist-3 operators. In section 
3, we present the actual calculation of the anomalous dimension matrix for the twist-3 
operators by employing standard QCD perturbation theory. The calculation is performed 
with the Feynman gauge. The loop integration is dimensionally regularized, and the minimal 
subtraction (MS) scheme is adopted. The details of the calculation will be discussed in the 
Appendices to make the discussion transparent. In section 4, we will discuss experimental 
implication of our result. 



2 Basic formulation 

2.1 Twist-three operators for hi 

In this section we shall briefly summarize general aspects of the chiral-odd spin-dependent 
distributions hi and hi relevant for our analysis. For the detail, we refer the readers to 
ref. 0. The QCD factorization theorem tells us that a cross section for an inclusive hard 
process can be decomposed into the perturbatively calculable hard cross section and the 
parton distribution function |T^. The latter is known to be written as the light-cone 
Fourier transform of the quark (or gluon) correlation function in a hadron. The chiral-odd 
parton distribution functions (renormalized at the scale n) in our interest are defined as 
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follows: 



+hL{x,ix^)M{p^ny -pun^){S ■ n) + hsix, fi^)M{S±f,n^ - S^yU^ 



{21) 



where \PS) is the nucleon (mass M) state with its momentum P and spin S {P^ = M^, 
S'^ = — M^, P-S = 0). We introduced the null vectors p and n by the relation P^ = p^ + ^n^, 
p"^ = = 0, p ■ n = 1, n'^ = p~ = 0, which specify the Lorentz frame of the system. The 



hght-cone gauge n ■ A = was employed in (|2.1| ). hi and hi are directly accessible by 
measuring the proper asymmetries in the polarized Drell-Yan process {h^ is twist-4 and 
is irrelevant for the following discussion.) Taylor expanding the bilocal operator appearing 
in the correlation function, l.h.s. of (|2.1|) , one can derive the relation between the moment 
of these parton distribution functions hi^i and the local operator, 

^M^w -Mn _ Sn^ii^cy^^iD^"' ■ ■ ■ iD^'-ijj, (2.2) 

where the covariant derivative = — igA^ restores the gauge invariance and the sym- 
bol Sn symmetrizes the indices /ii,...,/i„. Here and below, we often suppress the exphcit 
dependence on the renormalization scale of the local operators and the parton distribution 
functions. For the study of the twist-2 and -3 distributions hi and h^, it suffices to consider 
the piece of fft^'^t^^'^t^" which is further symmetrized among u, fii, ...^n- We thus introduce an 
arbitrary light-like vector (A^ = 0) and consider 

^/..A = 0A«-^--'^"A,A^, •■•A^„. (2.3) 

9^ ■ A can be decomposed into the traceless part 9 and the remainder T, 

0^A = 0>•A + T„^■A (2.4) 

by the condition that 

g.J^''^'-^- = guj^""^'-^- = g^,^,9^^''^'-^- = 0. (2.5) 

Qfj.ufj.i- fj.„ contains all the twist-2 effect in 9^"^^ "^'^ and it is related to the moments of hi{x, fi) 
as 



MnMfi)] = a„(/i), (2.6) 
{PS\9^- A{fi)\PS) = ('s'^p^+i - P^^P" + ^i-M^A^SP"-^) , (2.7) 

Ivj. \ Til I ^ / 

where we introduced the shorthand notation Ain[h{n)] = J dxx^h{x, fi) and k = k ■ A for an 
arbitrary four vector /c^. Also, we obtain for the moments of h^: 

Mnlhl] = -^Mn[hi] + Mn[hL]. (2.8) 



3 



The first term shows that /il receives a contribution from the twist-2 distribution hi. This 
piece is an analogue of the Wandzura-Wilczek contribution for g2. Mnlhi] of (W^ ) 
are directly related to the matrix elements of the twist-3 operator ■ A of ( p.4|) . Explicit 
calculation gives 



■ A 



A n-1 



+ 2 ^ „ 



(2.9) 



where d = iD ■ A. ■ A can be recast into the following form by using the relation 
[D^, Dy] = —igG^y with G^i, the gluon field strength tensor 



Ttt-A 



n 



n + 2 

where 0-A = 0^i-^"A 
1 



-A^Er, ■ A 



-^A^Nn -A-A^yll- 



Rnl ■ A, 



(2.10) 



Ml 



A^^ . Here the first term (the "EOM operator" ) is defined as 
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traces. 
(2.11) 

This operator vanishes by the naive use of the QCD equation of motion {ip — mq)ip = 0. 
We can set it to zero when we take its matrix element with respect to a physical state (such 
as the nucleon state) [|T6|, |1^], which is why it is discarded in ref.Q. However, this is not an 
operator identity and the mixing between En and the other twist-3 operators defined below 
should be taken into account during the course of renormalization. The second term is given 
by 

jyMi-Mn ^ Snmgifj'j^-f^'HD^'^ ■ ■■iD^'^i) - traces. 



Rn\ in the third term of ( p.lOj ) is defined as 



^J,l■■■^J,n 



n + 1 



(2.12) 



(2.13) 



D/il ■••Mil 



igC'^ ■ ■ ■ iD^"x/j - traces, 



(2.14) 



which explicitly involves the gluon field strength tensor, suggesting that the twist-3 operators 
truly represents the effect of quark-gluon correlations. By the combination of ^p" '^" in the 
form of the r.h.s. of ( |2.13| ), R'^i "^" can have definite charge conjugation property. 

With these definitions and the relation (|2.1CI|) , the n-th moment of the genuine twist-3 
piece of hi can be written down in terms of the nucleon matrix elements of the twist-3 
operators (see (p.8|)) 



Mn[hi 



n rrir 



n + 2 M 



Mn-i[gi] + Mn[hl]. 



(2.15) 



The first term is due to the contribution of the second term of eq.( ^.10D , and thus shows the 
quark mass effect; it is in fact the quark mass times the twist-2 operator corresponding to 
the (/i-distribution: 



(P5|5„7/^7^^752D^^ ■ ■ ■ iD^'-iPlPS) = 2Mn-i[gi{f^)]Sn{S^'' P^^ . . . _ traces). (2.16) 
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The second term designates the contribution from 

21 



Mn[hl] = E 1 



1=2 



n + 2 



with 



{PS\KT^''il^^)\PS) = 2bn,iif^')MSniS>''P''' . . . P^" - traces). 



(2.17) 



(2.18) 



By inverting the moments, ( |2.8| ) and ( p.l5| ) give the relation between the structure functions 
themselves: 



hL{x,fi) = 2x dy+-f 

Jx M 

for X > 0.Q 



X 



+ hl{x,fi^), (2.19) 



2.2 Renormalization of the twist-three operators 

We now proceed to discuss the renormalization of the twist-3 operators. As discussed in 



n+l 



N„ ■ A and E„ 



section 2.1, [{n + l)/2] + 1 twist-3 operators, Rn,i ■ A (/ = 2,3, 

A formally participate in the ra-th moment of hL{x,Q'^). It has been known that as 
a complete basis of higher twist operators one can always choose "canonical" ones which 
(1) are traceless and symmetric with respect to all the Lorentz indices and (2) have no 
contracted derivatives. Any noncanonical operators which could appear through radiative 
corrections can be transformed into the canonical ones modulo EOM operators by use of 
the relation [D^, D„] = —igG^^, and the physical matrix elements of the EOM operators 
vanish. These canonical operators mix with each other and also with the EOM operators 
under renormalization. Therefore we can choose Rn,i ■ A, Nn ■ A, and ■ A as a basis for 
renormalization. We further recall that the renormalization of composite operators generally 
involves the mixing with gauge noninvariant EOM operators which do not exist in the original 
basis [O]. We will come back to this point in the next section. 



The scale dependence of the physical matrix elements, e.g., &n,i(/^) of (|2.18| ) are determined 
by the anomalous dimensions of the corresponding composite operators. To see this, we 
write down the renormalization group equation for these operators. The bare- (Of) and the 
renormalized- {Oi) composite operators are related by the renormalization constant matrix 



(2.20) 



where Oi symbolically refer to i?^^"^", E^^'" '^", and A^^^'" '^". The renormalization group 
equation for Oi{ii) is obtained by using the fact that the unrenormalized operators do not 
depend on the renormalization scale: 



^hL{x) for a; < should be related to the "antiquark distribution" h^lx) as /il(x) 
conjugation. 



(2.21) 

-hLi^x) by charge 
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where the anomalous dimension matrix for {R,E,N} is defined as 

7.(/^)=/^^^^r.^(/^)- (2.22) 
In the leading logarithmic approximation, this equation is solved to give 

j 

where is the QCD running coupling constant, /?o = (H ~ (2/3)A^/) /(47r)^ and 70 are 

the lowest order (1-loop) coefficient of the /3-function and the anomalous dimension matrix 

7^J{9{^^))=lU^^r + 0{g{f^r). (2.24) 

In order to compute the Z^-factor of (|2.20|) , one imbeds the composite operators Oi 
into an appropriate Green functions with a convenient kinematics, and computes radiative 
corrections to this Green function. For a familiar case of the twist-2 operators, for example, 
the two-point functions with the on-shell external lines are usually considered. In the present 
case, however, the EOM operator E!^^'"'^^ should be retained as a nonzero quantity and the 
mixing of them with the other operators should be consistently taken into account: This 
means that it is necessary to compute Green functions with the off-shell kinematics for the 
external lines. 

Thus, in order to compute renormalization of Oi = for example, we may consider 

the truncated three-point Green function Fi{p, q,q — p) defined by 

Fiip, q, k){2nY5\p + k- q)G{p)G{q)D{k) 

= j d^xd^yd^ze'P''e-"iye'^'{T{0.ii){x)i){y)A,,{z)}), (2.25) 

where G and D are the quark and gluon propagators, respectively. (We suppressed the 
Lorentz and the spinor indices for simplicity.) We consider the three-point Green function 
for the off-shell quark and gluon external lines (not a physical state!), and therefore Fi with 
Oi = En do not vanish. In the next section, we present the one-loop calculation of this 
function to get Zij. 




0,(/i^ 



(2.23) 



3 Anomalous Dimension Matrix for Twist-3 Operators 

In this section we present the computation of the anomalous dimension matrix for the twist-3 
operators for h^. The calculation is performed up to the one-loop order. The loop integrals 
are dimensionally regularized and the MS scheme is employed for renormalization. Thus 
we keep only the simple dimensional pole proportional to 1/e in the one-loop amplitudes 
{e = {4 — d)/2 with d the space-time dimension) . We use the Feynman gauge for the gluon 
propagator, but the results should be independent of the gauge. 

As we discussed in the last section, we imbed the relevant twist-three composite opera- 
tors into the three-point function Fi{p,q,k) of ( |2.25|) assuming the off-shell kinematics. In 
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this case, the basic ingredients are the tree level vertices for the operators Rn,i, En, Nn 
corresponding to the diagrams shown in Fig. 1, which we call the "basic vertices". 
The three-point basic vertex of Rnj ■ A shown in Fig. 1 (a) becomes 



while those for En ■ A, Nn ■ A are given by 



(3.1) 



r(3) 



g_ 

2 



i=2 



n,ii 



m, 



i=2 

n 

,^7A^75AEp'"'^""'^^ 

i=2 



(3.2) 
(3.3) 



where k = q — p and t°- (a = 1, A'^^ — 1) is the color matrix normalized as Tr(t'^t'') = 

We first imbed Rn^i ■ A into the three-point function Fi of (p.25|) . The Feynman diagrams 
which give the one-loop radiative corrections to the operator Rn^i ■ A are the one-particle- 
irreducible diagrams shown in Fig. 2. In order to write down those amplitudes, the vertex 
for Rn^i ■ A corresponding to Fig. 1(b) is necessary in addition to the usual FejTiman rules: 



r 

2 



^/''^=tV°^AA^75A^w"-'g' 



n-l 



n-l+1 



i=2 



1-2 



+ Yl cr-'A^zj,t-tY-'A,ip + ky-''+'-%-kg^, + k^A,)r-' 

j=n—l+3 

+ {fi ^ u, k ^ k' , a ^ b)] — {I ^ n — I + 2) 
with k' = q — p — k. We have to add the counter term contribution (see ( |2.20| )), 



(3.4) 



z^,'Z2^/z,z, - 6u'j 7^S, ^ + zriZ2^z,z, {s^'l + {Zn, - 1) ATif)} + zr^z,^z,z,z„,Nil 

(3.5) 

to the sum of all the one-loop amphtudes of Fig. 2, and require that the total results be finite 
as e — s> 0. In ( |3.5| ), Z2 and Z^ are the usual wave function renormalization constants for the 
quark and gluon fields, while Zg and Zm are the coupling constant and mass renormalization 
constants defined by 



9 = ^f^-'g""; 
^9 



mn 



(3.6) 



where "B" denotes the unrenormalized quantities similarly to ( ^.2(J| ). Note that ^Z^n\^l^ ^n}^i 
and A/'P^ are proportional to the coupling constant and thus each term of ( |3.5| ) involves the 



factor Zg as a coefficient. Similarly, the factor Z^ appears because A/"^^^ contains the quark 
mass m. In the MS scheme, (the finite part of) the renormalization constants are chosen so 
that the counter term contributions ( p.Sj) precisely cancel out the terms proportional to the 
Xje pole from the dimensionally regularized one- loop amplitudes. 

For the case of the three-point functions imbedding ■ A or iV„ ■ A, we can proceed in 
a similar manner by interchanging the roles of Rn,i ■ A with ■ A or ■ A. These results, 
together with the well known results for Z2^3^g^rn (in the Feynman gauge), 

g"^ I g'^ g"^ 

^ ' -Cf\ ZgdZ'i = \— Cg] Zm = 1 — 3 Cf (3.7) 



(47r)% (477)%^' (47r)% 

completely determine the relevant Zj^-factor = 2,..., ,E,N). {Cp = ^^fWT ^"^^ 
Cg = Nc are the Casimir operators of the color gauge group SU(Ai'c).) 

The actual computation of the Feynman amplitudes of Fig. 2 is rather cumbersome due 
to the complicated structure of the vertices, e.g., (|3.1|) - (|3.4|) . Moreover, the participation 
of the gauge- noninvariant EOM operators would further complicate the computation (see 
below). In this light it is convenient to employ the following procedure (l)-(3) to make the 
computation tractable: 

(1) We introduce a vector ^2^, which satisfies the condition Cl = A^f2^ = 0. The contrac- 
tion of the implicit Lorentz index fi of Fi ( (p.25|) ) with fi^ kills off many terms and simplifies 
the computation enormously: For example, the basic vertices in ( |3.1| ) and ( p.2[ ) become 

7^i^J • n = -|fi,or"^A,z75(g - p) [p^-'t^ - t« (3.8) 

and 

£f ■ n = |^],a""AA^75 (j)""' + r^') t\ (3.9) 



This contraction brings another favorable effect: As was exemplified in ref. in the context 
of gT, the gauge-noninvariant EOM operators generally mix through renormalization in 
addition to the gauge- invariant EOM operator En [0 . Typically, those gauge-noninvariant 
operators are obtained by replacing some of the uncontracted covariant derivatives contained 
in i?^^ "'^" by the simple derivatives. Therefore, for large n, a large number of different 
gauge-noninvariant operators are expected to come into play. Due to this phenomenon, it is 
extremely difficult to identify the tensor structure of the one-loop amplitudes by the basic 
vertices ( ^.1D -( PT3D and those for the gauge-noninvariant EOM operators. On the other hand, 
after the contraction with Q^, we do not distinguish En- A and the gauge-noninvariant EOM 
operators; i.e., the contracted basic vertex ■ of ( |3.9| ) and the contracted basic vertices 
for the gauge-noninvariant EOM operators coincide due to the condition f2 = 0. It is easy 
to see that such an identification between the gauge-invariant and gauge-noninvariant EOM 
operators does not affect the prediction of the Q^-evolution of the moment sum rules, because 
of the property that the physical matrix elements of the EOM operators vanish. 

(2) For the computation of the three-point functions imbedding Rn^i ■ A and En ■ A, 
we set TJiq = 0. Taking this limit is legitimate for the present case employing the off-shell 
kinematics for the external lines, and amounts to neglecting the contribution of the basic 
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vertex A/"^^^ of ( p.3|) in the preceding discussions of this section. Clearly, this procedure still 
gives the correct results for the renormalization mixing between and En, i.e., for Zij 
with i,j = 2, ■ ■ ■ , [(n + l)/2], E. Actually, we need not compute the one-loop correction to 
the three-point functions imbedding En ■ A: The property {PS\En ■ A\PS) = immediately 
implies = (i = 2, • • • , [(n + l)/2], N). 

(3) In order to obtain the other components of Zij, it is sufficient to consider the two-point 
functions shown in Fig. 3 with the insertion of the relevant operators. We again employ the 
off-shell kinematics for the external quark lines, but the computation is now performed for 
the nonzero quark mass. The basic vertices corresponding to the operators En ■ A, Nn ■ A 
are given by 



n-l 



(3.10) 
(3.11) 



while the one for Rn.i ■ A vanishes. Thus, the computation of the one-loop corrections to 
the two-point functions imbedding Rn^i, En and Nn, by using the vertex given in (|3.1| )- 
( |3.3| ) and by following the steps similar to the case of the three-point functions, gives the 
components ZiE, Zi^ {i = 2, ■ ■ ■ ,[{n + l)/2], E, N). For the remaining components Z^i 
{i = 2, ■ ■ ■ ,[{n + l)/2]) we need not perform any actual computation: Those vanish because 
Nn is a twist-2 operator multiplied by a quark mass. 

We note that the several components of Zij can be obtained by different methods, giving 
a consistency check of our procedure: ZiE (^ = 2, ■ ■ ■ , [(n -|- l)/2]) by calculating the two- as 
well as three-point functions; Zen, Z^e by the calculation of the two-point functions and by 
a special property of the operators En, Nn discussed above. For all those cases, the different 
methods gave the identical results. 

In the above discussion, we completely neglected the flavor structure of the operator. 
Even for the case of flavor- singlet combinations of hi^i, there exist no gluon distributions 
which mix with them. In fact, if there would be any mixing between flavor- singlet hi and 
a gluon distribution, it would arise from the diagrams shown in Fig. 4. But all of them are 
identically zero because of chirality. The situation for hi is completely the same. 

Now we get all the Zjj-factor for the twist-3 operators. (The contributions of the relevant 
three- and two-point Feynman amplitudes corresponding to the diagrams of Figs. 2 and 3 are 
presented in Appendix A.) We summarize the final result in the following matrix form: 



If we express Z^j as 



ZimijJ-) ZiEifi) Zmifi) \ / -Rn,m(/i) 
ZEEifi) Enifi) 

Z^iv(/i)y V NM 



l,m = 2,- 



n + 1 



.(3.12) 
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then Xij is given as follows: 



n + 1 



,E,N 



(3.13) 



X, 



Im 



m + 1 



G 



\n 



I + 1]2 [/ - 1]^ 
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+2 



+ 



\n — l + 2 Ij n — l + 2 — ni l — m 



I (21-1 _ m(/ - 3)\ _ 1 /2n - 2/ + 3 _ m{n - I - lY 
^2 [[l-l]2 ~ [I - 1]3 J ~ 2 [[n-l + l]2 ~ [n-l + l]3 ^ 



+{Cg-2Cf) 
+2(-l) 



-| \n—l—m (A \l—m 

--LJ n-l^m-2 \~^) 1-2^1- 



n-l + 2-m) n-m+lCl-l {I - m) n-m+lQ-r 
l-\Cm-l n-l+lCm-l 



[Z-l]3 [n-Z + l]3, 



(2 < m < / - 1), 



Cr 



■2(n-2/ + 2)[n-/ + l]2 " 2 ^"^"-'^^ + ^'"^ + ^""'^^ + 



1 f{2l-n-2){n-l-l) 1 + 2 1-3 \ 1 
^2 I, [n - Z + 1]3 ^ Z(Z-l) ^ [Z-lJaJ ~ T 



+2(Cg - 2Cf) 
+(-1) 



-1)' 



+ 



-1)' 



[n-l + 1]3 



/-I 



2(n-/ + l)(n-2/ + 2) 

— Ci? (25';_i + 2Sn-l+l — 3) , 



Im 



G 



1 1 n - 2m + 2 

n — I + 2 — m l — m 2[n — / + 1]2 

{n - I - l){2m - n - 2)' 



-(Cg - 2Cf) 



2[n-l + 1]3 

2( — 1)™ {n-l+lCm-l — ("l)" n-;+lC'm-z) 



n—l^m—l 



[n-l + IJa 



n—l—m rt 



{m - I) m-lCm-l {n - I + 2 - m) n-m+lCl-l 

"n + 1 

Z + 1 <m < 



XiE — 2Cf 



,[l]2 [n-l + 2]2j' 



X, 



IN 



-4Cp 



}l-l]s [n-l + l]3j' 
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XeE — 2(1 — Sn)CF, 



(3.19) 



X 



NN 



n(n + 1) 



n+l 



(3.20) 

E,N), 



where ^„ = E"=i 1/j and [j]^ = j(j + l) ■ ■ ■ U+k-l). With these X,, (2, j = 2, 3, ^ 
the anomalous dimension matrix for the twist-3 operators Rn,i, En and take the form of 
the upper triangular matrix as 



lij 



(3.21) 



From (|2.23| ), the Q^-evolution of the nucleon matrix elements are given by 

-X/167r2/3o' 



[^1 

E 

m=2 



C^n(Q') 



)2N \ -^/167r2/3o 
^ \ -Xjvjv/167r2/3o 



Im 



^a(/x2 



rf„(/i2)(3.22) 



IN 



(3.23) 



,/M)A^„_i[(7i(/i2)] (see (|2.16| )). This is our main result of this 



where we set dni^fi"^) = (m, 
paper. 

Before leaving this section, we shall briefly discuss the Q^-evolution of the twist-2 dis- 
tribution hi{x,Q^) because it appears in as a Wandzura-Wilczek analogue ( iW^ ) and 
( 2.19| )). There is only one operator 6!^ ■ A (eg. ( p.4D ) for each n in the twist-2 level, and thus 
there is no complication arising from operator mixing. The lowest order coefficient of the 
anomalous dimension can be obtained from the one-loop diagram of the two-point functions 
shown in Fig. 3. We note that the contribution from the diagram in Fig. 3 (a) vanishes, and 
the contribution to the Z-factor from Figs. 3 (b) and (c) is the same as those for fi and gi 
from the same diagrams. By using a similar technique as above, we obtain for the Z-factor 
for the composite operator ■ A, which is now a single constant and not a matrix: 



9' 



-CpK 



n+l; 



(3.24) 



with 



j=2J 



By substituting the result into (|2.22|) , we get the anomalous dimension 

^2 



7 



9 



n+l- 



This governs the Q^-evolution of the first term on the r.h.s. of ( |2.8| ) (see eqs. ([27 
( |2.7| )). The same result was obtained in ref.0| from the Altarelli-Parisi equation for hi. 



(3.25) 

(3.26) 

and 
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4 Examples of Q^-Evolution 



Here we present some examples of the Q^-evolution of the twist-3 distribution hL{x, Q'^) by 
using the result obtained in the previous section. Since there is no mixing of with the 
gluonic distribution, we shall consider a distribution for one quark flavor. We also neglect 
the contribution from the quark mass operator in the following. 

For the 3-rd and 4-th moments oi Hl, only one twist-3 operator ; contributes. Their 
anomalous dimensions are (ignoring the common factor g^STT^) 104/9 and 1099/90 for the 
3-rd and 4-th moments, respectively. This gives for A^3,4[/iL((3^)] as 




These curves normalized at ^ = 1 GeV are shown in Fig. 5. Here and below we set 
Nf = 3 and Aqcd = 0.5 GeV. For comparison, we also plotted the moments of the twist-2 
distributions fi and hi. 

From this figure, one can clearly see that the third moment of the twist-3 distribution evolves 
significantly faster than that of the twist-2 structure function. 
For n = 5, the anomalous dimension becomes the 2x2 matrix: 

/ _ 202 191 \ 

^=[71 _i49 , (4.3) 
^ 60 180 ^ 

and the eigenvalues of 7 = -X/{11 - |iV/) are 1.435 and 2.005. We thus get for A15[^l(<5^)] 
as 



'a{Q 



2\ \ 1-435 




MnihLiQ')] = (0.41665,2(/i) + 0.19365,3(/^)) ^ 

+ (0.01365,2(//)- 0.05065,3(//))(^J . (4.4) 

In principle, if one measures Ai^lhLiQ^)] at two different values of with sufficient accuracy, 
one could fix the two matrix elements &5,2(a*) and b^^sdi), and the measurement of A15[/il(Q^)] 
at different supplies a test of the QCD evolution. Since we do not have any physical insight 
on these matrix elements, we plotted M-slhiiQ'^)] normalized at // = 1 GeV in Fig. 6 with 
four moderate values of A(/i) = 65,3(/i)/&5,2(/^) = —4.0, —2.0, 1.0,4.0 at /i — 1 GeV. One can 
see from the figure that there is a large variation in the Q^-evolution among different choices 
of X{n). This fact suggests that a nucleon model and a nonperturbative technique of QCD 
can be tested by comparison of their prediction on A(/x) with future experiments. 

As a measure of the large Q^-behavior of the moments, we have plotted in Fig. 7 the 
lowest eigenvalues of the matrix 7;^ (l,m = 2,.., [^^]) with Nf — 3 for hi, Qt (defined 
as the twist-3 part of gx), hi, and /i as a function of n. For qt, we used the result for 7 
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obtained in ref. [T^. (The authors of ref. [0 calculated the anomalous dimension matrix of 
twist-3 operators only for the even moments of qt, since they discussed it in the context of 
the deep inelastic scattering.) From this figure we expect that the moment of the twist-3 
distributions evolves faster than that of the twist-2 distributions. This fact indicates that 
the measurement of qt and hi greatly serves as a new test of QCD at high energies. If 
one looks into more detail of fig. 7 one sees that the chiral-odd structure function, hi and 
hi, evolves slightly faster than the chiral-even ones with the same twist, /i, gi and qt- In 
general, however, actual form of the Q^-evolution in a finite window strongly depends on 
the relative magnitude among bn,i{[j) as we saw in the above example for n = 5. We thus 
should take the result shown in Fig. 7 only as a rough measure for the asymptotic behavior 
of the Q^-evolution. 



5 Summary and Conclusion 

In this paper we have studied the Q^-evolution of the moments of the chiral-odd twist-3 
spin structure function hL{x,Q'^) in the standard QCD perturbation theory. For the n-th 
moment of the twist-3 part of /^^^(x, Q^), M-riyhhiQ'^)], independent twist-3 operators, 

Rn^i (/ = 2, ■ ■ -, ), Nn and En play roles. The Q^-evolution of M-nlhiiQ"^)] is governed by 
the anomalous dimension matrix for the operators iV„, We thus have calculated 

the one-loop correction to the three-point Green function which imbeds these operators. 
Although the physical (on-shell) matrix element of the EOM (equation of motion) operator 
En vanishes. En mixes with Rn^i and Nn through renormalization and it is essential to take 
into account this operator mixing to determine the anomalous dimension matrix. In order to 
incorporate the mixing correctly, we employed the off-shell kinematics for the external lines. 
Using the minimal subtraction (MS) scheme we obtained the renormalization constants in 
the one-loop order in the form of the upper triangular matrix which really shows the mixing 
among {Rn^uNn.En}. As an example of the Q^-evolution of hi, we have studied the 3-rd, 
4-th and 5-th moments of hL{x,Q^) in detail for rriq = 0, comparing them with the known 
twist-2 distributions, /i, gi and hi, and the other twist-3 distribution qt- To consider the 
case rriq = is a sufficiently good approximation for the u and cZ-quark distributions, and 
is a reasonable one for the s-quark distribution, since the nucleon matrix element of iV„ 
is expected to be small for the s-quark. The notable features of the Q^-evolution of these 
moments can be summarized as follows: 



1. The Q^-evolution of the 3-rd and 4-th moments can be predicted uniquely, since there 
is only one relevant operator Rn^2 {n = 3,4). (Note the nucleon matrix element of 
En vanishes.) Compared with the twist-2 distributions /i, gi and hi, the moments of 
hL{x,Q^) evolves significantly faster. 

2. The fifth or higher moments of hi receives contribution from two or more operators, 
and thus the Q^-evolution depends on the ratio among these matrix elements at a 
reference scale fi. They can, in principle, be determined by measuring hL{x,Q'^) at 
several values of Q^, and the measurement at different values of gives a test for the 
QCD evolution. 
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3. The lowest eigenvalues of the anomalous dimension matrix for each moment of /il and 
gx (the twist- 3 part of qt) are much larger than those of twist- 2 distributions, which 
implies that the Q^-evolution of these twist-3 distributions is significantly faster than 
the twist-2 distributions. These numbers for are slightly larger than those for gx, 
which suggests that the chiral-odd distribution hx evolves faster than the chiral-even 
one gT- The similar tendency has been known for the twist-2 distributions, hi and gi. 

The anomalous dimension matrix for the twist-3 operators obtained in this paper de- 
termines only the Q^-cvolution of the moments of hL{x,Q'^). In order to predict the Q^- 
evolution of the whole x-dependent distribution hx^x, Q"^), we need to construct a generalized 
Altarelli-Parisi equation for a relevant multi-parton distribution function. This is because 
a higher-twist distribution essentially represents a correlated quark-gluon distribution and 
Q^) is only a particular projection of this generalized multi-parton distribution. This 
work is under way and will be published in a future publication. Nevertheless, we already 
found in this work some peculiar features in the Q^-evolution of the moments of hi,{x,Q'^), 
which we hope will be measured in the future collider experiments. 
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Appendices 



A. One- loop Feynman amplitudes 

In this Appendix we present the expressions for the one-loop Feynman amphtudes for the 
truncated one-particle-irreducible Green functions with the insertion of the composite oper- 
ators Oi = Rn,i ■ A (/ = 2, ■ ■ ■ [(ra + l)/2]), ■ A and N^, ■ A. 

First we consider the three-point function with the insertion of Rn^i ■ A. We set rriq = 
following the discussion of (2) in sect. ^. The Feynman diagrams for the one-loop corrections 
are shown in Fig. 2. The Feynman amplitudes are contracted by and are grouped to 
possess a definite charge conjugation. (Note that, in the Feynman diagrams of Fig. 2, (d), 
(f), and (h) are connected by the charge conjugation to (c), (e), and (g), respectively.) Each 
group can be expanded by the basic vertices TZ^l^ and (Here and in the following, we 
denote TZl^^n ' ^ ^.nd Sj^j^^ ■ Q of ( P^ and ( pi?D simply by TZ^n}n Sj^j^K A sample calculation 
of the diagrams (e) and (f) are given in Appendix B.) 

Fig. 2 (a) gives 0. 

Fig. 2(b) gives: 



9' 



-C, 



G 



E 



(/ + m)(m — 1) {n — I + 2 + m){m — 1) 



\2(/ -m)[l -1]2 2{n-l + 2-m)[n-l + l^^ 

(n + 2)(/- 1) 



[{n+i)/2] /(^2n-l-m + A){n-m+l) 



Si — Sn-l+2 + 777 + o7 I , n\ 

21 2[n — 1 + 2) 



+ E 

m=l+l 



2{m - l)[n - I + 1]2 



I 7? (3) 

2(n-2/ + 2)[n-/ + l]2j 

(n — / + 2 + m) (m — 1) 
2{n - I + 2 - m)[n - I + 1]2^ 



Fig.2(c) + Fig.2(d) gives: 



9' 



{2Cf - Cg) 



E 2(-ir f- 

m=2 \ 



0/ ( n-l+lCl-l - (-1 



n~l+lCm-l l-lCm-1 

|n.-/ + l]3" |;-1]3, 



[(n+l)/2] 



2f-ll 



^ [n-i + ih 

( 1 1 



7^(3) 



\[l]2 [n-l + 2]2^ 
Fig.2(e) + Fig.2(f) gives: 



(n-/+lCm-l ~ (~1)" n-l+lCm-l) T^n^m 



(A.2) 



9' 



■ i-i 

E 

.m=2 



{2Cf - Cg) 



-1 



\n—l—rn 



n-lCm-'2 



(—1) l-2Cl-rn 
{I — m) n-m+lQ-m (n — / + 2 — m) n-m+lQ-l, 
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G 



I n-l + 2' ^ ' ■'"•"^ 



- {2Cf - a 



7^(3) 

-mi-i] 



(n - 2/ + 2)(n - / + T) + ^^^^^'"^ + ^""'"^^ " + C'cy | 7^S 



m=l+l 

(y " ;y^TT2 

Fig. 2(g) + Fig.2(h) gives: 



m—l r~< 



1 Xn— i— m 



^(m-l) ^^iCrn-i {n-l + 2- m)„_„+iQ_i 

r(3)" 



7^(3) 



(A.3) 



G 



[n 



1 /(/-3)(/ + l-m) / + 2 



.S n - 1]3 



/-l)(n-/ + 3 



n-l + A 



\n 



1-3 



l + lh 
I + 2 



1 

+ 2 V[/-l]3 ' [l-lh 



[n 



[n-l + l]2/ ■■ 
/- l)(n-2/ + 2) 



(n-l-l){n-2m + 2 



- E 

m=l+l 
1 



/ + 1 



2[n-/ + l]3 

^ ^ ^-(3) 



- / + 1]3 
7^(3) 



/ + 3 



(A.4) 



The coefficients in these expansions diverge as e = (4 — (i)/2 — > with d the space-time 
dimension. By adding the counter term contribution ( p.5[ ) to the sum of eqs.( |A.lD -( [A.4| ) and 
by requiring that those counter term contributions cancel the 1/e pole terms, we obtain the 
renormalization constants given by in the MS scheme. 

Next we proceed to calculate the two-point function with the insertion of Rn^i - A, En - A 
and Nn ■ A, following the discussion of (3) in sect. |[ We keep the quark mass nig as a 
nonzero quantity. The computation can be performed in a similar manner as in the case 
of the three-point functions. The three-point vertices necessary to compute the one-loop 
diagrams of Fig. 3 are given by (|3.1|) - (|3.3| ). Again the amplitudes combined appropriately 
can be expanded by the basic vertices ( p.lO| ) and 

The one-loop correction to the two-point function with i?„ ; ■ A comes from Figs. 3 (b) 
and (c). It gives 



9' 



\[l]2 [n-l + 2]2^ 
4 4 ' 



V[/-l]3 [n-l + l\,^ 
For the one- loop correction to the one with En ■ A, Fig. 3 (a) gives 

IGvr^e n 



(A.5) 



(A.6) 
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and Fig.3(b)+Fig.3(c) gives 



9 



1 




Finally, for the one-loop correction to the one with ■ A, Fig. 3(a) gives 

.2 2 



9 



-C, 



167r^e n{n + 1)' 
while each of Figs. 3 (b) and (c) gives the same contribution, 



9' 



n 1 \ 
J=2 J J 



(A.7) 



(A.8) 



(A.9) 



By adding the appropriate counter term contributions, we can determine the counter 
terms in the MS scheme, which give the results ( |3.17| )-( |3720| ). Note that ( p.l7| ) are obtained 
from the three-point functions as well as from the two-point functions, giving a consistency 
check of our methods. Also, Zen = can be verified explicitly by using the results ( A. 6 ) 
and ( |A.7| ) (see the discussion of (3) in sect. H). 



B. Sample calculation 

In this appendix we describe the details of the calculation of the one-loop Feynman am- 
plitudes. We choose the diagram (e) and (f) of Fig. 2 with the insertion of Rn,i ■ A as an 
example; the other diagrams can be calculated using the similar technique. 

First we write down the Feynman amplitude JF(e) for the diagram (e) by using the vertex 
and the usual Feynman rule (in the Feynman gauge): 



abc^c 



[{iP - kT'^(i'^ -iP- A;)'"'g""'} if 

' n-l+l l-l 

E (p - ky-Y'^^'-^t' - E(p - ky~'p^-'-^r-' } (<? - P)t''t^ 



+ { E {p - k)"^\q - ky-^-^'^^q''-^ 

j=n-l+3 



n 



E - k)^~^iq - ky-^-^q""-^ ) {q - py't^ 

3=1+1 



1^' 



(B.i; 



where q and p are the incoming and the outgoing off-shell quark momenta, and k is the 
gluon loop momentum. We work in the massless quark limit {rriq = 0) following (2) of sect. 
3. The Lorentz index /i corresponding to the external gluon line is contracted by fi^, which 
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kills off the terms involving in the vertex ( p.4|) . Also, owing to this property, the factors 
involving the gamma matrices can be easily evaluated: 



(B.2) 



We follow the standard procedure: We use the Feynman parameterization to collect all the 
denominators of the quark and the gluon propagators. After shifting the /c-integration, many 
terms can be dropped by the condition = 0. We retain only the divergent contribution 
of the ^-integration: 



IbiT^e Jo 



(n-l+1 l-l 
j=2 j=2 



+ c 



r-'\{q-p)). (B.3) 



where t = (1 — x)p, s = —xp + q. We perform the Feynman parameter integral by 



dxi^-H^-^-^ 



1 i-l-' 



r=0 



l+rCr 



-P 



(B.4) 



The resulting formula of JF^g) involves double summation. These terms can be simplified 
into single summation by interchanging the order of summation and by using the relation 
T,r=m rCm = n+iCm+1- Then, wc obtaiu 



G 



(-1) 

m=2 ~ ^) n-rn+lCl 



1 



167r%2 ^ ^ [ 2 \V ' n-l + 2 

+Cf {{S^-i+i - l)p--'t' - (Si^i - l)p'''r-'} (q-p) 

2 

n f -| \m—l 

m=l+l ~ ^) m-lCm-l 



pn-l+l-l-2 



(B.5) 



The amplitude corresponding to the diagram (f) can be computed in a similar manner 
as above. The result should be related to ( [B.5| ) by charge conjugation. By adding this result 
to ( p.5|) , we obtain as the total result: 



(/) 



X 



2^V-^^uO-^^AAi75t'' 
167r%2 ^ 



n-l + 2 



^pn^l+l^l-2 j^pl-2qn-l+l 



+Cf(S,_, + S„_,+i - 2) (p"-'q'- - p'-'r ) (4 - P) 



-1) 



l—rn 



" (-1)— 



-m=2 ~ n-m+lCl- 



pn~m^m—2 pm—2^n—rn 



(q-p) 



(B.6) 



This should be expressed as a hnear combination of the basic vertices ( p.8| ) and ( p.9| ), which 
we denote simply by TZ^^] and This is readily performed for the terms proportional to 
Cp or Cp — Cg/'^- For the terms proportional to Cq, we use the identity 



p 



p 



m=2 



It 



(B.7) 



m=i+l 



Now we obtain 



-ic.-^)i't 



2) K-;; 



(3) . 



Ca (I 



m=2 ~ n-m+lC*Z-m 

2U-/ + 2 



7+1 ~ m-lC*! 



m=l 
7^(3) 



m—l 



(B.8) 



Here 7^(f)„ with m = [(n + l)/2] + 1, ■ ■ ■ , n appear. These can be expressed by those for 
m = 2, ■ ■ ■ , [(n + l)/2] using 'R-^n^n-m+2 = ~'TZn}m,i we obtain ([A.3|) of Appendix A. 
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Figure Captions 

Fig. 1 (a) Three-point basic vertex for Rn,i, En and A^„. (b) Four-point basic vertex for 
Rn^i necessary for tlie calculation of the diagrams shown in Fig. 2. 

Fig. 2 One-particle-irreducible diagrams for the one-loop correction to the three-point Green 
function {eq.{^n^). 

Fig. 3 One- loop corrections to the two-point function relevant for the calculation of Zie, 
ZiN. These diagrams are also used for the calculation of the anomalous dimension of 
the twist-2 distributions. 

Fig. 4 Diagrams which could cause mixing between a fiavor-singlet quark distribution hi 
and a gluon distribution (if any). These diagrams are identically zero for the chiral-odd 
distribution. 

Fig. 5 The Q^-evolution of the 3-rd and 4-th moments of /^^(x, Q^) normalized at /i = 1 
GeV. The 3-rd moments of twist-2 distributions /i and hi are also plotted for compar- 
ison. 

Fig. 6 The Q^-evolution of the 5-th moment of Q^) normalized at /z = 1 GeV for four 
moderate values of A(/i) = —4.0, —2.0, 1.0,4.0. 

Fig. 7 The smallest eigenvalues of 7 = — X/167r^/5o as a function of the dimension of the 
moment, n. 
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